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Abstract. The simultaneous effect of both disorder and crystal-lattice pinning on the equilibrium 
behavior of oriented elastic objects is studied using scaling arguments and a functional renormaliza- 
tion group technique. Our analysis applies to elastic manifolds, e.g., interfaces, as well as to periodic 
elastic media, e.g., charge-density waves or flux-line lattices. The competition between both pinning 
mechanisms leads to a continuous, disorder driven roughening transition between a flat state where 
the mean relative displacement saturates on large scales and a rough state with diverging relative 
displacement. The transition can be approached by changing the impurity concentration or, indi- 
rectly, by tuning the temperature since the pinning strengths of the random and crystal potential 
have in general a different temperature dependence. For D dimensional elastic manifolds interact- 
ing with either random-field or random-bond disorder a transition exists for 2 < D < 4, and the 
critical exponents are obtained to lowest order in e = 4 — D. At the transition, the manifolds show a 
superuniversal logarithmic roughness. Dipolar interactions render lattice effects relevant also in the 
physical case of D — 2. For periodic elastic media, a roughening transition exists only if the ratio 
p of the periodicities of the medium and the crystal lattice exceeds the critical value p c = 
For p < p c the medium is always flat. Critical exponents are calculated in a double expansion in 
/! = p 2 /Pc — 1 and e = 4 — D and fulfill the scaling relations of random field models. 

PACS. 68.35.Ct Interface structure and roughness - 71.45.Lr Charge-density-wave systems - 
64-70. Rh Commensurate-incommensurate transitions - 68.35.Rh Phase transitions and critical 
phenomena - 05.20. -y Statistical mechanics - 74-60. Ge Flux pinning, flux creep, and flux-line 
lattice dynamics 



£^ ' 1 Introduction the magnetic field in the case of Ising magnets, respec- 



tively. 

The thermal roughening transition (RT) from a faceted a third closely related example for the RTs are the 

to a smooth surface of a crystal is one of the paradigms lock-in transitions of periodically modulated structures 

of condensed matter physics ||. It has been observed di- into an underlying crystalline matrix. Examples for the 

rectly for a number of materials, in particular on surfaces modulated structures are spin or charge density waves 

of crystalline Helium 4 @. A similar transition occurs be- (CDWs) in anisotropic metals §§, mass density waves 

tween a rough delocalized and a flat localized state of the in sup erionic conductors §, polarization density waves in 

interface in the three dimensional Ising and related lattice incommensurate ferrolectrics § , Wigner crystals @,mag- 

models g@. At the RT temperature T R the free energy of netic bubb i e arrays @, flux line lattices (FLLs) @|2) 

a step on the surface vanishes. Apart from its shape also an d nem atic elastomers fl|, etc. The phenomenology of 

other physical properties are influenced by the presence of tnese bulk systems i s even more rich than that of sur- 

the RT, e.g., the growth (drift) velocity u of the crystal or interfaces, since in general a large number of different 

surface (interface) changes dramatically at the RT from locked-in structures separated by incommensurate phases 

u ~ exp(-C//) for T < T R to u ~ / for T > T R f2|. may occur . 

Here / denotes the driving force density, which is propor- .... 

, . ij.il j -xr a r,i i -i . ; • i r? Common feature of the examples mentioned is the ex- 

tional to the difference An of the chemical potentials of . . ^ . .. 

the crystal and its melt in the case of crystal growth and f tence ° f a ^-dimensional deform-able object (i.e. D = 1 

tor single flux lines, D = 2 for surfaces or interfaces, 

a Email address: te@thp.uni-koeln.de D = 3 for CDWs or the FLL etc.) the energy of which 
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- in the absence of lattice effects - can be described by 
elasticity theory. The distortion of the elastic object is ex- 
pressed in general by an TV-component displacement field 
4 (x),i = 1...N, where N = 1 for surfaces, interfaces and 
CDWs, N — 2 for flux lines or FLLs etc. These elastic 
objects interact with the potential of the underlying crys- 
talline matrix, which favors (in a commensurate situation) 
certain spatially homogeneous ("flat") configurations be- 
cause of their lower energy E. On the other hand, thermal 
fluctuations favor a spatially inhomogeneous ("rough") 
configuration because these have the higher entropy S. 
The real macroscopic configuration is that of the lowest 
free energy F — E — TS, from which we expect the occur- 
rence of a roughening transition at a temperature T = Tr. 

Later we will subdivide these elastic objects into mani- 
folds, which comprise single flux lines or domain walls, and 
periodic media, i.e. periodic lattices of flux lines or domain 
walls, CDWs etc. This distinction will become necessary 
because manifolds and periodic media couple in different 
ways to the disorder. However, a distinction is unnecessary 
as long as we consider merely thermal fluctuations. 

The interest in the investigation of the roughening 
transition goes also beyond its immediate applications in 
condensed matter physics. In computer simulations of con- 
densed matter systems a fake periodic lattice potential - 
in addition to that of the crystalline matrix - appears fre- 
quently in the form of a grid on which the simulations 
are done, which has nothing to do with the physics of 
the problem. In order to retrieve physically relevant infor- 
mations from these simulations one has to make sure to 
be above the roughening transition of this fake potential. 
This remark applies also to lattice gauge theories in high 
energy physics [ 14|,|l5| . 

Besides of its frequent occurrence the second interest- 
ing aspect of the thermal RT is its universality. Although 
the RT temperature Tr itself is non-universal and de- 
pends in particular on the clastic properties of the object 
under consideration (e.g. on the orientation of a vicinal 
surface relative to the main crystalline axises), the uni- 
versal features like the singular behavior of the thermo- 
dynamic quantities are the same for all models with the 
same spatial dimension D and number of field components 
N. Moreover, the thermal RT for D — 2-dimensional sur- 
faces was shown to be in the same universality class as 
the metal-insulator transition of a 2D-Coulomb gas uvL 
the Kosterlitz-Thouless transition of a 2D XY- model Jl7[ 
(both with an inverted temperature axis), the phase tran- 
sition of the 2D sine-Gordon model jl8| and the transition 
of the ID XXZ-chain at T=0 (T|. 

D = 2 is indeed the upper critical dimension for the 
thermal RT since D > 2 dimensional manifolds are flat 
even without a lattice potential. 

Interesting but less known is the fact that Kosterlitz 
pijf and Forgacs et al. |2l[ considered the thermal RT for 
1 < D < 2 interface dimensions and found a superuni- 
versal (i.e. dimension independent) logarithmic roughness 
at the transition. Finally, the thermal RT disappears in 
D = 1 dimensions: 1-dimensional interfaces are rough at 
all finite T 



Disorder is an even more efficient source of roughen- 
ing of manifolds than thermal fluctuations and physically 
important as well j23fl . 

Roughening of interfaces due to disorder was consid- 
ered to determine the lower critical dimension of the ran- 
dom field Ising model |p4| , p5|j2l|| and the mobility of do- 
main walls in disordered magnets |^ , ^8|j29| , |30| , |3l| , ^2| . It 
was argued that in the presence of disorder interfaces of 
dimensions D < 2 are always rough, but undergo a rough- 
ening transition as a function of the disorder strength for 
D > 2 dimensions ]33|,|4|]. Bouchaud and Georges [|5) 
considered explicitly the competition between lattice and 
impurity pinning in 2 < D < 4 dimensions using a vari- 
ational calculation. Surprisingly they found three phases: 
a weakly disordered flat, a rough glassy and an interme- 
diate flat phase with strong glassy behavior. For D < 2 
only the glassy rough phase was expected to survive, but 
that cross-over regions of the two other phases would still 
be seen in the short scale behavior. However, it should be 
kept in mind that the variational calculation is an uncon- 
trolled approximation and typically reproduces the results 
obtained from Flory-like arguments. 

It is the aim of the present paper to give a unified 
renormalization group description of the roughening tran- 
sition of elastic manifolds and periodic media in the pres- 
ence of quenched disorder. Some preliminary results were 
already published |3^|37]]. In this paper we focus on the 
derivation of the renormalization group flow equations and 
their solutions for the different cases. To make the paper 
more self-contained some of the results reported already in 
Refs. [^6 37 are included here as well. The physics is com- 
plicated by the occurrence of various input length scales 
like the periodicity of the crystalline matrix oq /p, the cor- 
relation length of the disorder £o and the periodicity of 
the periodic medium I. 

The outline of the paper is as follows. In Sec. 2 the 
model for the pinned elastic objects is described. The dif- 
ferent classes of disorder correlations are briefly reviewed. 
Section 3 analyses the model on the basis of simple scaling 
arguments. A qualitative picture of the phase diagram is 
developed and the characteristic length scales are identi- 
fied. In Section 4 we proceed with the derivation of the 
functional RG equations to lowest order in D = 4 — e 
dimensions. The functional RG flow of the disorder corre- 
lator is reduced to a flow of single parameters in Section 5. 
The fixed points of this RG flow are analyzed for vanish- 
ing lattice potential, reproducing the known results. The 
effect of a finite lattice potential is considered in Section 
6. The critical exponents of the resulting roughening tran- 
sition are calculated to lowest order in e for the different 
types of disorder. Experimental implications of our results 
and related recent experiments are discussed in Section 7. 
The last section concludes with a summary and discus- 
sion of our results and briefly comments on related open 
questions not treated in this paper. 
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2 The Model 

As mentioned already in the Introduction, the elastic ob- 
jects are described by an ./V-component phase field 0(x) 
with components l (x), i — 1, . . . , N, which interacts both 
with an underlying periodic crystalline lattice potential 
Vl{4>) and a random potential Vr(0, x) specified below. 
The total Hamiltonian is then given by 

H = 1 J ^x{~V0- V0-V L (0(x))-V JI (0(x),x) 

(1) 

The strength of the pinning potentials is measured in 
units of the elastic stiffness constant 7 which controls 
isotropic elastic deformations. More complicated elastic 
moduli could be considered but turn out to be unessen- 
tial as long as we are interested in universal properties. 
For manifold models, the elastic object is embedded into 
a larger space of dimension d = D + N whereas for peri- 
odic media the dimension of the space in which the elastic 
structure is embedded is equal to the internal dimension, 
d = D. The different realizations of the model ([l]) are 
summarized in Table [l| 

The periodic potential Vl(0) can be rewritten as 



y L (0) = ^W). 



(2) 



i=i 



if the original crystalline matrix is assumed to have the 
symmetry of an TV-dimensional hyper-cubic lattice with 
spacing oq. Below we will measure all length scales in units 
of ao, i.e. we will put ao = 1. The resulting V{4> 1 ) is then a 
periodic function of periodicity 2w/p. The random pinning 
potential Vr(0, x) is assumed to be Gaussian distributed 
with a vanishing average Vr(0, x) = and two-point cor- 
relation function 



Va(0, x) W, x') = R(</> - <j>')5^ (x - x'). (3) 

The starting point for our further calculations is the 
replica Hamiltonian 7i„, which follows, e.g., from the con- 
sideration of the disorder averaged free energy 



F = -T lim - (Tre 
n 



-1 



(4) 



with 
7~C n 



E 

a/3=l 



d U X 



2 

_L 

2T 



V0 Q • V<j> - V L {<f> a ) 

R(<J> a - /3 ) 



0a(3 



(5) 



The actual form of the disorder correlator i?(0) depends 
on the type of disorder and on the model under consider- 
ation. 

We give here a brief discussion of the main classes of 
disorder correlators: 

For elastic manifolds one has to distinguish between 
random field and random bond systems. 



(i) For random field system 
correlator is given by 



Ro{<t>) 



f 1016 for |0| »£> 
n ° \ 2 for |0| « Co 



4U25H the bare disorder 



(G) 



Here \fSo denotes the strength of the random field, Aq 
is proportional to the concentration n; mp of the random 
field impurities. £0 is of the order of the maximum of the 
intrinsic domain wall width and the correlation length of 
the random field. The result (Q) was derived in p4p^ ] for 
N = 1 component fields. Its formal generalization to arbi- 
trary N > 1 is straightforward but without direct physical 
application. However, this case will be included here in or- 
der to enable the comparison of our results with that of 
the variational approach of Bouchaud and Georges [ p5| 
which is expected to become exact for N — > 00. 

As an other interesting realization of a sine-Gordon 
model with random field disorder, we suggest surfaces of 
crystalline 4 He in an aerogel. Here the silica strands act 
as quenched impurities on the crystal surface. Fractal cor- 
relations in the aerogel density appearing on short scales 
should be negligible if the correlation length £|| of the sur- 
face exceeds all structural length scales of the aerogel. 

(ii) For random bond systems with short range correla- 
tions in the disorder the bare disorder correlator is given 

by 

i?o(0) = A>£o 2 exp(-0 2 /2£o)- (7) 

The length £0 has the same meaning here as for random 
field systems. Note that in our notation £g includes a fac- 
tor N such that £q/N remains finite in the limit N — > 00. 

(iii) In the case of periodic elastic media the form of the 
bare disorder correlator is dominated by the periodicity 
of the elastic medium. To show this more explicitly we 
start, for simplicity, with a scalar field <j). The periodically 
modulated density can be expanded in a Fourier series 



E 



p m cos[m(qx + 0(x))] 



(8) 



where q denotes the modulation vector of the locked-in 
phase. The coupling to the disorder potential U(x) is of 
the form 

7 R (0,x)=C/(x)p(x,0) (9) 
where U (x) = Uq 5(x — Xj) — n lmp ). Here x, and 



j imp 



are the positions and the concentration of the impurities. 
The disorder average is performed by integration over all 
possible impurity positions, 



.4 



n 



d u x, 
V 



A({ Xl }). 



(10) 



Here the integration extends over the volume V of the sys- 
tem. With (B), (B), (0) and neglecting rapidly oscillating 
terms we get then in a straightforward manner 



M<t>) 



1 



00 



P^O 2^ Pm COS ( m ^)- 

=0 



(11) 
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Table 1. Different physical realizations of the model (|l|); topological defects are excluded. 



Elastic object 



d 



D N 



Interface |i||||§6| 



density wave 

Xy-magnet 
in crystal and 
random field 

flux line lattice [hlLM 



4 He surface 
in aerogel 

nematic elastomer 
in external and 
random field |l3f 



> 2 d- 1 1 



2-Ku/a 
u displacement 



> 2 d 1 Sp(x) = Pi cos (qx + 0) 



> 2 d 1 



> 2 d 2 



> 2 d- 1 1 



> 2 d 1 



arct an (S'y/S'z 



2nu/l 
I flux line spacing 

2nh/a 
h height variable 

&rct&n(n y / n x ) 
n unit director field 



q = Q/p 



> 1, number 
of easy directions 



p = //a 



— 10| (random field) 
e -0 /S (random bond) 



Eq#o #Q cos(Q0) 
Q reciprocal lattice vectors 



This expression can be extended to vector fields </>(x) as 
discussed by Giamarchi and Le Doussal Q in the context 
of flux line arrays in superconductors, and gives 

R (cp)= ]Ti? Q cos(Q0) (12) 

Q#0 

with Fourier coefficients Rq where Q are the re- 
ciprocal lattice vectors of the periodic medium. Gener- 
ally, rewriting the periodic density p(x, <p) in terms of 
a displacement field <fi, also gradient terms of the form 
— po Yli=i ^i^ 1 are generated in Eq. (||) which lead to terms 
of the form di^djcjP^ in Eq. (||). But these gradient terms 
can be shown to be strongly irrelevant for the large dis- 
tance behavior for D > 2 |39|] . Therefore, in the following 
this interaction will be neglected. 

3 Scaling arguments 

In this section we present some elementary scaling consid- 
erations for the system described by the Hamiltonian (|l|) , 
which give a qualitative picture of the phase diagram, the 
critical dimensions and the characteristic length scales. 
Since we are interested only in qualitative estimates, we 
neglect all numerical factors of order unity. 

We start with the elastic manifolds and restrict our 
analysis to the interface case N = 1. The extension to 
N > 1 does not change the qualitative conclusions derived 
for N = 1. Since disorder fluctuations turn out to be much 
stronger than those of thermal fluctuations, we put T = 
from the very beginning, postponing the discussion of the 
influence of finite temperatures to the end of this section. 

In the absence of disorder, Ao — 0, the interface is 
flat. Typical excitations consist of terraces of height 2n/p 




Fig. 1. Two dimensional manifold forming a terrace to gain 
energy from the random potential without any bulk energy 
cost from the periodic potential. 



and linear extension L. The terraces are bound by steps 
of width w ~ l/py/v and energy -E s t op — I^/vqL ' 1 /p 
where we used for simplicity Vl{4>) — vq cos pcj) for the lat- 
tice potential, see Fig. [l]. Disorder can lower the energy 
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of terraces since the locally displaced interface sees a dif- 
ferent random potential. To be specific, we find for the 
energy of a hump of total height <j> consisting of 4>p/2tt 
steps 



Etotaih ' 



(13) 



The second part represents the energy gain from the ran- 
dom field. Minimizing the total energy ( |l3| ) with respect 
to </> we obtain pip ~ (L/Lr) 2 ~ d with 



L, 



(v /£ A p) 



1/(2-25) 



(14) 



for D < 2, and L has to be replaced by the lattice spac- 
ing ao for D > 2. In deriving (14) we implicitly assumed 
£o ^ a o/p with ao = 1 for the estimate of the disorder 
energy in jl^). In the opposite limit £o ^ a o/p pertur- 
bation energy gives for the disorder energy \J A L D . A 
useful interpolation formula for Lr is therefore 



vq(£oP+ 1) 
A £oP 



l/(2-D) 



(15) 



For D < 2 the flat phase is unstable and the interface is 
therefore rough on scales L > Lr even in the presence 
of the crystal potential. Since the energy cost for terraces 
vanishes now on length scales L ~ Lr, the system is de- 
scribed on larger scales by an effective elastic Hamiltonian 
with a stiffness constant j c g. A rough estimate for 7 e g fol- 
lows from balancing the elastic energy and the bare energy 
cost for a terrace of height 2tt/p on the length scale Lr, 



JefiL R 



D-2 




p, 



(16) 



which gives 7 c ff ~ jL R /w |3j^|34|]. On scales L L R the 
hump height scales as pip ~ 2tt\L/Lr) i > where £ = (4 
D) /3 is the roughness exponent for random fields |24|,^5| . 
Lr is physically meaningful only if it is larger than the 
bare kink width, Lr > w. This condition can be rewritten 
as 

vo > v 0tC « {-j^) • (17) 

If «o ~> vq, c +, 7 e ff changes into 7. For v Q < Vq iC , the lattice 
potential can be neglected and vq has to be replaced by 
«o,c in the above formulas. 

So far we have considered the case D < 2. For D = 2 
the length scale Lr can be obtained from an argument 
which has been originally applied by Morgenstern et al. 
p(| and Binder |^6| to the surface of an Imry-Ma domain 
in the two dimensional random field Ising model. Consider 
again a terrace of height 2?t/p and size L in the otherwise 
planar interface. As can be seen from Eq. (|l^), both the 
step energy and the random energy scale proportional to 
L. To calculate Lr one has to consider the distribution 
of random energies -Edis of such a terraces resulting from 
different disorder configurations. The mean value of this 
random energies vanishes, but its variance increases with 
L. For large enough terraces, the distribution for Edis be- 
comes Gaussian, 



V(E d 



1 



AnL^Ao^o/p 



exp 



1 



2 8tt£MoWp 



(18) 



Since the total energy cost for a terrace of circular shape 
is given by the sum of the elastic energy E c \ = 2ir 2 ^/vqL / p 
and -Etfis) the probability V < that the total energy becomes 
negative can be easily estimated to be 

V<= f 61 dE dis V{E dis ). (19) 

J — OO 

This integral evaluated in the limit A — » gives 



1 



P^qA q 



■ exp 



4 p£oA> 



(20) 



Hence the density V < ~ (Lr)~ 2 of such terraces tends 
exponentially to zero as A — > 0. Neglecting algebraic fac- 
tors, terraces of arbitrary size appear spontaneously be- 
yond the exponentially large length scale 



j R ~ exp 



vo 



p£oA 



(21) 



such that the interface becomes asymptotically rough. Thus, 
the flat phase is unstable to a weak random field in D < 2 
dimensions. 

For D > 2, we have to expect from the previous con- 
siderations that the interface becomes rough only for vq < 
Uo,o i- e -i that there is a RT as a function of vq (or Aq) at 
vo ~ vo, c . Since one could object that the argument for 
the occurrence of a critical value vq jC for vq is restricted to 
D < 2, we give here an additional argument for the exis- 
tence of a roughening transition for D > 2 by comparing 
the pinning forces which originate from the periodic and 
the impurity potential, respectively. If there is no disorder, 
an external force (density) has to overcome a threshold 
value 

fv w Pivo (22) 

to depin the interface. On the other hand, for i>o = the 
threshold force from the impurity potential on a length 
scale L is given by Q 



(23) 



where we introduced the Larkin-Fukuyama-Lee length 
La ~ (Co/^o) 1 ^ 4 ^ and the roughening exponent ( of 
the manifold in the absence of lattice pinning. The tran- 
sition between the lattice pinned flat and the impurity 
pinned rough phase occurs if both force densities are equal 
on the length scale of the step width w. This gives 



v,^p- 2 {p^{A /a? /{i - D) . 



(24) 



For random field systems, C = (4 — -D)/3 if p£o <C 1 
(asymptotic regime) and C = (4 — D)/2 if p£o 1 (per- 
turbative regime), in complete agreement with (|l7|). 

For completeness we additionally show the stability of 
the rough phase to a weak crystal potential. Indeed, lowest 
order perturbation theory in i>o yields an effective Hamil- 
tonian where the periodic crystal potential is replaced by 
a mass term m 2 2 /2. Assuming a Gaussian distributed 
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field 0(x), which is justified to lowest order in e = 4 
we obtain on length scales L > La for the mass 



D, 



in = p vocos(p4>) = p «ocxp[— p </> 2 /2] 

= p 2 v exp[-27T 2 (L/L A ) 2C ]. (25) 

Comparing this exponential behavior in L with the L 2< »~ 2 
scaling of the elastic energy, we conclude that the per- 
turbation is always irrelevant since C > for the elastic 
manifold models. Likewise, the flat phase is stable to a 
weak random potential as can be shown easily by per- 
turbation theory with respect to Vr. From the stability 
of both phases in D > 2 dimensions we can expect that 
there is a roughening transition for D > 2 as a function 
of v a (or Aq) at v ~ v 0tC (A ). 

In the rough phase, the mean square displacement of 
the interface can be estimated again from the energy 

"total 



cf + ^L^^^-^A^L^' 2 . (26) 



Here we have included both the elastic and the step energy 
term to describe the crossover from vq <C «o,c to vq — v c . 
Moreover, we have treated the lattice pinning term in a 
self consistent harmonic approximation to account for a 
situation in which p<f) > 2n. In the rough phase, the hump 
height asymptotically scales with the roughness exponent 
( = (4 - D)/3 as a result of the balance between the 
first and the last term in Eq. (|26|). However, there is an 
intermediate length scale region up to £|| in which the 
roughness of the interface is dominated by the competition 
between the second and the last term of Eq. ( [26]) . In this 
region the roughness increases only logarithmically, 



62 „ ^ {D _ 2)ln f _L 
P 2 \Lr 



(27) 



This result will be confirmed below by the RG calculation 
which yields also the exponent v\\ of the crossover length 
£|| to first order in e = 4 — D. Finally, for vq 3> «o,c the 
interface is always flat. 

So far we have treated only the case of random fields. 
For random bond systems similar arguments apply but the 
asymptotic roughness cannot be obtained correctly from 
a Flory argument. 

For periodic elastic media, arguments identical to those 
used before for manifolds give a lower critical dimension 
D c {= d c ) = 2 for the occurrence of a roughening transition 
and the estimate ( |l7| ) for the critical value of vq jC if we set 
SoP» I- 

To estimate the influence of a weak crystal potential 
in the rough phase, we compare again the scaling of the 
elastic energy with that of the effective mass m 2 as given 
by Eq. (|25|). In what follows, the logarithmic roughness of 
the structure, <j) 2 = ^(4 — D) \u(L/La), on length scales 
L > La and for 2 < D < 4, cf. Sec. [f| is crucial. Due to 
this weak roughness, the effective mass decays only alge- 
braically on large length scales L, 



p v L 



-V 

L A ) 



The elastic energy simply scales as L 2 and thus we con- 
clude that the crystal potential is a relevant perturbation 
if 

P<P ^^t=^ (29) 

leading always to a flat phase. For p > p c the rough phase 
is again stable and one has to expect a roughening tran- 
sition. Possible subtleties about exchanging the limits of 
large L and small e have been neglected here, but this 
result will be reproduced below by a detailed RG double 
expansion in e = 4 — D and fi = p 2 /p 2 — 1. 

So far we have neglected thermal fluctuations. In D > 
2, where we have anticipated the occurrence of a rough- 
ening transition, thermal fluctuations are irrelevant. In 
D < 2 fluctuations are still dominated by disorder fluc- 
tuations on large scales. However, for weak disorder and 
sufficiently high temperatures T > Tr, a length cross-over 
from thermal to disorder dominated pinning phenomena 
will occur as briefly discussed in p3,p4| . 



4 Functional Renormalization Group 

Next we employ a functional RG calculation, which starts 
with the disorder averaged replica-Hamiltonian (|^). We 
apply a momentum shell RG to this Hamiltonian, in which 
the high-momentum modes of 0(x) in a shell of infinites- 
imal width, |q| £ [Ae~ dl ,A], are integrated out. The mo- 
mentum cutoff A will be kept fixed by rescaling coordi- 
nates and fields according to 



e <v, 



0(x) = e^V(x') 



(30) 



Since the behavior of the elastic objects should be de- 
scribed by zero-temperature fixed points, we allow the 
temperature to be renormalized by keeping the elastic 
stiffness constant 7 fixed. From the above scale changes 
we obtain the RG flow of the temperature, 



^ = (2-D-2C)T. 



(31) 



(28) 



It can be shown that this flow equation is exact for van- 
ishing lattice potential due to the tilt symmetry of the 
Hamiltonian (|) for V L {<j>) = ||. As we will see below, 
a finite lattice potential leads to a renormalization of the 
elastic stiffness which can be rewritten as a renormaliza- 
tion of the temperature. This leads to an additional term 
in Eq. ( |3l| ) but does not change the temperature into a 
relevant RG variable. 

The RG equations for Vl{4>) and R{4>) will be derived 
in their functional form following the method developed by 
Balents and Fisher jl3) . The two last terms of the replica- 
Hamiltonian (||) will be considered as a small perturba- 
tion 7i P ort of the elastic gradient term which we denote 
by TIq. To eliminate the high- momentum modes, the field 
is divided into slowly and rapidly varying contributions, 
</>(x) = </> < (x) + > (x). The feedback from integrating 
out </> > (x) can be calculated by expanding the effective 
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replica-Hamiltonian for the slow fields (x), 

H n = - In J [d0>] cxp(-H n ), (32) 

in powers of 7i pC rt- The first terms of this cumulant ex- 
pansion are given by 

T~tn = Ho + (W pe rt)>,c — 2(^pert)>,c + g(^ P crt)>,c — 

(33) 

Here (. . .)>, c denotes the connected expectation value over 
> (x) with respect to Tto. Generally, new interactions 
which are not present in the original Hamiltonian are gen- 
erated by the RG transformation and have to be included 
in Tipcrt if they are relevant in the RG sense. We postpone 
the discussion of such sort of terms to the step in the RG 
calculation where they appear for the first time. 

Now we start to evaluate the series (H|) order by order. 
To obtain the feedback of first order from the perturba- 
tion, we rewrite H pcr t in its Fourier representation, 



H 



pert 



i W e *^(x)£ (k) 

A1 a p Jx > k 



e-^-^Vb(k), 



(34) 



x.k 



where we use here and in the following the vector k for 
Fourier transforms with respect to (/>, whereas the vector 
q will be used to denote momenta in x direction. Fur- 
thermore, we have introduced the abbreviations J , = 

fd D xfd N k/(2n) N , Q/3 (x) = a (x)-^(x) and 7 /T - 
1/T. 

After decomposition into slow and rapid fields, the av- 
erage over the fast modes in Eq. ( |34| ) can be easily per- 
formed using Wick's theorem. For simplicity, we will drop 
the label < on the slow fields in the averaged expressions, 



(Wpert)>,c 

-y 

2T 2 ^ 



a/3 



e *k0 Q(3 (x) e -k 2 TG>(O)(l-5 Q0 )^ k ) 



+ -E / e 4k ^W e -^ k2TG >(°)V L (k). 
T a Jx.k 

Here we have introduced the free two-point function 

A 

c 

2 



(35) 



G>(x)= / -j= J (T^p^T^ (36) 

\q\=Ae~ dl 



where the label > on the integral denotes integration over 
the infinitesimal momentum shell. In the limit dl —> this 
function, at x = 0, becomes 

G> (0) = (2tt)- d S D A°- 2 dl = K D A°- 2 dl, (37) 

where So is the surface area of a unit sphere in D dimen- 
sions. Since we have to calculate the expectation value of 



the perturbation to first order in dl, the exponentials in 
Eq. ( |35| ) depending on the two-point function G>(0) can 
be expanded with respect to dl. This yields 

— (7"^pcrt)>,c 



a/3 



K D A D ~ 2 dl 

2T 



didiR{(f> a p) -ndARiO) 



a/3 



(38) 



The corresponding diagrams are shown in Fig. Here 
and below the partial derivatives di act on the field <f>, 
and repeated Latin indices are summed over from 1 to 
N. Whereas the first two terms are simply the original 
perturbation, the remaining new parts are generated by 
the RG. The first and the last term of this new part is 
of the appropriate form to renormalize R(4>) and Vl{4>) 
but they are reduced by a factor of T compared to the 
original ones. Thus they are irrelevant at the zero temper- 
ature fixed points. The second new term proportional to n 
contributes only to the renormalization of the free energy. 
Therefore we conclude that in first order of the cumulant 
expansion there is no renormalization of R(4>) or Vl (</>)■ 



Fig. 2. Diagrams arising in first order of the cumulant expan- 
sion. The closed lines represent G>(0). Both diagrams con- 
tribute to the free energy only, see Eq. (pq). 



Now we turn to the analysis of the second order contri- 
butions to the renormalizcd Hamiltonian. Averaging over 
the fields d> > , this contribution becomes 



1 



9 C^pcrt) >,c 



/ / 1 V P ,:kl ^« 1 ^( Xl )+ lk2 ^c 2 3 2 ( X2 ) 

^ 2 /3 2 



xe -k 2 TG>(0)(l-,5 Ql3l )-k 2 TG > (0)(l-5 c , 2 3 2 ) 
xe -kikaTG>(xi-X2)(tfaia 2 +«jS 2 /3 2 -*c«ijSa-*/s 1 a2)_R(l Cl )^(l C2 ) 
j L V e ' k i^ ni /3 1 ( x i)+ jk2 ^ Q2 ( X2 ) e -k?TG>(0)(l-A-„ l/3 ) 

xe -ik 2 TG>(0) e -k 1 k 2 TG > (x 1 -x 2 )(5„ 1Q2 -5^ 2 )^ ( - ki ^ i( - k2 > ) 
j L_ V »k 1 cvi (x 1 )+ 4 k 2 </> Q2 (x 2 ) i(k 2 +k 2 )TG>(0) 

2T 2 ^ 

xe -k 1 k 2 TG > (x 1 -x 2 )5 Q1Q2 y i ( ki ^ i ( k2 - ) | _ d t j ^39) 
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where d.t. stands for disconnected terms which have to 
be subtracted from the above expression. To obtain the 
contributions to first order in dl, we have to expand the 
exponentials depending on G> up to second order in T . 
Sorting out the disconnected parts, the terms proportional 
to G> (0) are canceled. Performing the summation over 
the replica indices and transforming back to real space, 
we obtain the following second order result, 

2 C^pert)>,c 

+^E 9 ^(<M x ))^(<Ux')) 

a/3 

E W0 q (xO)}g>(x - x') 

a ' 
^ x ' x ' ^ a/3 7 



a/3 



2d i d j R{4> af} {*))d i d j R(Qj\ 
l-Y, didjRWdidjVL (0 Q (x')) 



2T 



■^a i a i y i (0 Q (x))5 i 9 J -y i (0 a (x')))G 2 > (x - x' 



(40) 

corresponding to the diagrams in Fig. ^. The three terms 
within the first curly bracket have a contribution only at 
large momenta since G> (x) is composed only of momenta 
within the shell. The first of these contributions is a 3- 
replica term which does not fit into the form of the orig- 
inal replica-Hamiltonian (|^). Such terms with more than 
two replicas result from non-Gaussian correlations in the 
disorder and are generated by the RG. Anticipating that 
R(4>) will be of order e at the fixed points, the renormal- 
ization of R{4>) and Vl{4>) resulting from 3-replica parts 
will be of higher order in e since the 3-replica terms are 
generated only at large momenta. Similarly, one can ar- 
gue that also the other two large momenta terms in the 
first curly bracket do not influence the RG results to low- 
est order in the small parameters of the RG expansion. 
Note that it may become necessary to introduce a second 
small parameter which controls the fixed point value of 
the lattice potential Vl (</>)■ Therefore, we will drop in the 
following all terms generated by the RG which contribute 
only at large momenta. 

From the terms in the second curly brackets, the first 
one which has three replica indices would renormalizc terms 
with three replicas which would appear in the original 
Hamiltonian with a temperature dependence ~ T~ 3 such 
that the above term can be neglected at a T = fixed 
point. The next two terms are relevant and lead to a renor- 






Fig. 3. Diagrams arising in second order of the cumulant ex- 
pansion, where the internal lines represent G>(x — x'). The 
graphs with one internal line contribute only at large momenta 
and can be neglected. 

malization of R((f>) and V/, (</>), respectively. The last term 
is also of the right form to feed back into Vz, (</>), but it is 
down by a factor of T and hence negligible. 

To evaluate the two relevant contributions of Eq. ([if)]) , 
we note that main contribution to the x'-integration comes 
from x' = x since G>(x — x') then becomes maximal. 
Therefore, we make of change of coordinates according to 
x' — > x + x' and expand the interactions with respect to 
x' . This gives formally expressions of the form 

f ir(x)F(x')G? > (x-x')= f F(x)F(x + x , )G^(x') 

J x.x' J x,x ; 

= f F 2 ( X ) [ G 2 > (x')+ / F(x)VF(x). / x'G^x') 

■/ x J x' •/ x J x' 

+^(x)a x „9 x „f(x)£xXG 2 > (x')+0(x' 3 ), (41) 

where n and m are summed over from 1 to D. The part 
which is linear in x' vanishes due to the symmetry of 
G>(x'). Since G>(x) is composed only of momenta of the 
shell, the kernel G>(x) can be easily integrated to first 
order in dl, 



f G 2 >(x) 

J X 



= K D A D ~ i dl. 



(42) 



Generally, one may expect also the appearance of quadratic 
gradient terms from second derivatives of -F(x) in the 
above expansion, resulting in a renormalization of the elas- 
tic stiffness constant 7. Whereas the absence of any Vl 
terms in the relevant contributions ~ R 2 renders gradient 
terms strongly irrelevant, they have to be taken into ac- 
count for the RVl part. Up to second order in V</>, the 
non-vanishing parts of a gradient expansion of the penul- 
timate term of Eq. ( [40| ) are 

-o^E/ W(o){a f a,-VL(0 a (x)) 

+ ~d i d j d k d l V L (<l> a (x))x'Vct> k a (xi) ■ x'V0L(x)}G| (x'). 

(43) 

The first term contributes to the renormalization of Vl 
whereas from the last part one may expect a change of 
the elastic stiffness. Thus to lowest order in the gradient 
expansion, we end up with the following contributions, 

S^R = ^K D A D - i dl {d l d j R(4>)d l d J R{(p) 
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-2^9^(0)^9^(0)} 
S (2) V L = —KDA^dldidjRWdidjVLty), (44) 

where S^R and S^Vl denote here and in the following 
the n th cumulant order feedback to the renormalization of 
R(4>) and Vl (</>), respectively. 

To analyze the last term in Eq. (^3|) we note that 
didjVL(4>) ~ §ij since we have assumed the particularly 
simple form (pfi for the lattice potential. Therefore the 
diagonal structure of the elastic kernel is preserved. To 
obtain the relevant contribution to the stiffness constant 
7, the gradient term has to be expanded in terms of eigen- 
functions of the linear operator of the RG flow of Vl{4>) 
JiH . This linear operator is given by 

CV L = (2 - 2QV L + C<P%V L - ^didjR^didjVL, (45) 

where the last term stems from the second order contri- 
bution (]44|) and the remaining parts arise from rescal- 
ing. The determination of eigenfunctions is simple since 
didjR(0) = —SijA in all cases of interest. It is simple to 
see that the eigenfunctions are given by 



E m (4>) = COS^TOi^ 1 ) 



COs(pTOjV<?!> 



(46) 



with m = (mi,...,mjv) a vector of integers. Here we 
have taken into account the rescaling of </>(x) by allowing 
for a periodicity a of the lattice potential which depends 
on the flow parameter I such that da/dl = — £a. Then the 
eigenfunctions multiplied by the simple exponential factor 
exp(A m Z) are solutions of the linear RG flow dV^jdl = 
CVl- The scaling dimension of this eigenfunctions at a 
fixed point is given by the corresponding eigenvalues 



2C m 2 lim p 2 A. 

2 /— ^OG 



(47) 



It will turn out below that at these fixed points, where the 
lattice potential Vl(4>) has a finite value, the eigenvalues 
become A m = 2(1 — m 2 ) up to corrections of the order of 
the small parameters of the RG expansion, e.g. e. There- 
fore, we see that higher harmonics of the lattice potential 
with m 2 > 1 are strongly irrelevant for the large scale be- 
havior of the elastic object. This justifies the assumption 
of Eq. (J2j) for the lattice potential. 

To analyze the last term of Eq. (^), this term has to 
be written as a linear combination of the operators 



(48) 



For m 2 > 0, all these operators are strongly irrelevant 
with scaling dimensions A m — 2+2£ at the fixed points with 
finite lattice potential. Therefore only the projection of the 
gradient term in Eq. ( [43| ) onto the constant eigenfunction 
E (<p) can potentially contribute to the renormalization of 
the stiffness constant. But this projection vanishes since 
Vl{4>) is simply a sum of first harmonics. Thus 7 is not 



renormalized to first order in the lattice potential, and we 
have no contribution to 7 from second cumulant order, 



5W7 = 



(49) 



where 5^ denotes the n th order feedback to 7. 

To obtain a finite fixed point for Vl{4>), we have to 
calculate the RG flow up to the first non-vanishing non- 
linear order in Vl (</>)■ Therefore, we have to consider now 
contributions from third order of 7i pC rt- In this order only 
terms of the order RV? are relevant as will be explained 
now. The terms ~ R 3 generate only three-replica parts 
which will be of order e 2 and can therefore be neglected. 
On the other hand, the contributions ~ V^ which would 
be able to renormalize Vl are reduced by a factor of T 
and thus are also negligible at the zero-temperature fixed 
point. For the terms of the order R 2 Vl the discussion is 
slightly more complicated. The corresponding two-replica 
terms generated in this order are of the form 



1 

T 2 



J2 aidjR&^didkRfr^djdkVLfa), (50) 

a/3 



which does not fit into the form of the terms present in the 
original Hamiltonian. To analyze this unusual term, one 
has to expand again the lattice potential Vl(4>) into eigen- 
functions of the linear operator of the RG flow equation 
for Vl(4>)- This will enable us to determine the relevance 
of contributions of the form (|50|). For that purpose we 
rewrite ( p0| ) as a linear combination of operators defined 

by 

S ^ = ^i I J2 d ^ R ^)d l d k R(ct> aP )E m (cf, a ). (51) 

The scaling dimension of these operators can easily be ob- 
tained in terms of A m . It turns out that the operators S m 
for all m are strongly irrelevant with scaling dimensions 
-L>-2-2C + A m = -4 - 2m 2 + 0(e). Therefore, terms of 
the form ([30|) are generated by the RG but are negligible 
for the large-scale behavior. 

The only terms in third order remaining to be dis- 
cussed are proportional to RVj 2 . In contrast to the other 
combinations of R and Vl treated above, these terms pro- 
duce a relevant feedback to renormalize Vl- The corre- 
sponding part of (Hp Crt )>,c is given by an expression sim- 
ilar to those appearing in Eq. ( |39| ) . In third order the rel- 
evant contributions come from the T 3 terms of the expan- 
sion of the exponentials containing the function G>(x). 
Integrating out the rapid modes and transforming back to 
real space, we end up with the result 

/*! E^^ i? ( )^^ yi (^( X2 ))^ afe ^(^( X 3)) 

xG>( Xl - x 2 )G > (x 1 - x 3 )G>(x 2 - x 3 ) (52) 

represented by the graph in Fig. ||(a). We obtain in lowest 
order of the gradient expansion in Eq. (41) the following 
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third order contribution to Vl{4>), 

6 {3) V l = -^K D A D ~ 6 dld i djR(0)d i d k V L (<l>)d j d k V L (4>). 

(53) 

But for this sort of feedback also higher order terms of 
the gradient expansion are relevant. The corresponding 
expression in second order in V</> is given by 

2? /-iE ^dMo)aid k diV L (4> a {^i))d j a k a m v L {4> a ^i)) 

xx 2 V^( Xl ) • x 3 VC(xi)G > (x 2 )G > (x 3 )G>(x 2 - x 3 ). 

(54) 

As before, the diagonal structure of the elastic part is pre- 
served by this feedback due to didjVL ~ &y. In contrast to 
the second part of Eq. (|43| ) , this contribution is quadratic 
in Vl- Therefore, also the first harmonics of the lattice po- 
tential yield a non- vanishing projection onto the constant 
eigenfunction Eq{4>) since sin 2 (p</>) = (1 — cos(2p</)))/2. 
Performing the spatial integration over x 2 and X3 to low- 
est order in dl and e, we obtain the relevant feedback to 
the stiffness constant from third cumulant order, 

£ (3) 7 = -K D A D - s dld l d ] R(0) 

x VoididkdiVLi^djdkdtVL^)]. (55) 

Here Vq denotes the projection operator onto the eigen- 
function Eo((f>). 



(a) 



(b) 





Fig. 4. Relevant graphs which appear in third (a) and fourth 
order (b) of the cumulant expansion. In both cases the disorder 
vertex contributes only a ^-independent factor didjR(0). 



Below we prefer to keep the stiffness constant fixed by 
allowing the temperature to flow, so we have to rewrite the 
feedback (|55|) as a renormalization of T. But this means 
that also the feedback to the functions R{(f>) and Vl(0) 
is changed due to the overall factor 1/T in front of the 
Hamiltonian. The effective but non-rescaled quantities are 
then given to first order in dl by 

T cff = (1 - 5 (3) 7 )T, (56) 
R cS = (1-2S^ 1 )R + S^R, (57) 
VL.eff = (1 - ^ ] l)V L + 5 (2) V L + 6 {3} V L . (58) 

The part 7 is quadratic in Vl, and it is tempting 
to neglect this term in Eq. (|5^ ) since it produces cubic 



Vl terms which appear irrelevant compared to the term 
S^Vl which is quadratic in Vl- But it turns out that this 
quadratic term does not contribute to the RG flow of the 
first and only relevant harmonics of the lattice potential. 
The reason for this is that a squared first harmonic pro- 
duces only a constant and a second order harmonic after 
projecting onto the eigenfunctions E m (<p). 

As explained above, the lowest non-linear Vl term in 
the flow equation of Vl will be of third order. Since all rel- 
evant contributions have to include at least one disorder 
vertex ~ R, we have to consider also possible feedback to 
Vl from the fourth order cumulant. Simple counting of the 
R and Vl orders suggests that in this order terms propor- 
tional to R 2 V L 2 and RV 3 have to be taken into account. 
For the first type of feedback similar arguments as those 
for the contribution of Eq. ( 50 ) show that it is irrelevant 
in the RG sense. In contrast, the second kind of terms gen- 
erates a relevant renormalization of Vl ■ This feedback can 
again be calculated along the lines described above. Now 
the relevant connected low-momenta parts come from the 
T 4 terms of the expansion of the exponentials with respect 
to the two-point function G>(x). To first order in dl the 
4 th order contribution becomes 

5 {i) V L = -^KDA^dW^R^dkVL^djd^L^) 



x d k d m V L (cj}). 



(59) 



This part follows from the lowest order of the gradient 
expansion in Eq. ( |4l| ) ; higher gradient terms are irrelevant 
since they could lead to a renormalization of 7 only of 
third order in Vl and, therefore, could produce higher 
order terms only in the flow equation for Vl- Therefore 
contributions to 7 of higher than second order in Vl can 
be neglected. 

Moreover, it is important to note that in higher than 
fourth order of the cumulant expansion no relevant terms 
are generated which renormalize R and Vl to lowest non- 
linear order. For example, in fifth order terms appear which 
are cubic in Vl and of the right form to renormalize Vl 
but with a coefficient proportional to R 2 . Therefore we 
can stop our analysis in fourth order of the cumulant ex- 
pansion. 

Keeping track of all relevant feedbacks calculated above 
and introducing the function 

w(4) = didjR^didkdM^djdkdmVL^), (60) 

we obtain the resulting functional RG equations, valid to 
lowest order in e = 4 — D, 

^ = (2-D-2( + V {W(<j>)])T=-8T, (61) 



dR(4>) 
dl 



(4 - £ - 4C)i?(0) + C<^i?(0) 



-didjRWdidjRty) - didjRWdidjRM) 



dV L {(t>) 
dl 



2T {W(ct>)}R(ct>), 

(2-2C)Vl(0)+CW4(0) 



(62) 
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- -didjRiO) [didjVLW + d i d k V L (cl>)d j d k V L (<t>) 

+ d t d k V L {<t>)djd m V L {<t>)d k d m V L (</>)] 

+ V [W(<I>)]V L (<I>), (63) 

where we made the replacements A~ 2 Vl{4>) ~^ Vl{4>), 
K D A D - 4 R(ct>) -» i?(0) with K D X = 2 D - 1 n D / 2 r(D/2). 



5 Absence of lattice pinning: The rough 
phase 

Much is known about the equilibrium behavior of clastic 
manifolds and periodic elastic media embedded in a dis- 
ordered background for D just below the upper critical 
dimension four. In the following, we will briefly summa- 
rize the known results for the roughness exponent £ in the 
various cases. In addition, we will develop a simplified cal- 
culation scheme to obtain the effective disorder strength 
or, equivalently, the RG flow of the disorder potential on 
sufficiently large but finite length scales. This technique 
becomes useful if we study the competition between ran- 
dom and crystal pinning in Section || 

In the absence of lattice pinning, the RG flow of the 
disorder correlator R((f>) has been studied repeatedly in 
the past JIIIIS The RG flow is determined by Eq. 
(H) with V L {4>) = 0. The fixed point functions R*{4>) 
are always of the order e, but their functional forms may 
be quite different depending on the bare functions Rq{4>) 
at / = 0. There exist mainly three different locally stable 
fixed points, each with its own basin of attraction. The 
corresponding solutions of Eq. ( |62|) with Vl{4>) = will 
be analyzed in more detail for the different types of struc- 
ture and disorder in the following subsections. Within this 
analysis the main focus is on the simplified scheme which 
allows for a reduction of the functional flow equation to 
a flow of at most two parameters. This novel representa- 
tion will be helpful below for the inclusion of an additional 
crystal potential in the RG scheme. 



5.1 Elastic Manifolds 

For the elastic manifold models, we have to distinguish 
between long-range correlated and short-range correlated 
random potentials. Both types of correlations can be phys- 
ically realized for domain-walls (N = 1) in, for example, 
disordered ferroelectrics or Ising magnets. Whereas a ran- 
dom field leads to long-range correlated randomness, dilu- 
tion forms short range correlated random bond disorder. 
In the following, we will generalize our considerations to 
structures with arbitrary N . For the rest of this section 
we set Vl{4>) = which corresponds to W(<j>) = in Eq. 



5.1.1 Random Fields 

It is easy to show that the R(4>) ~ \4>\ behavior for large 
|0| of the bare disorder correlator is preserved under the 



RG flow. Therefore we have to search for rotationally in- 
variant solutions of the flow Eq. (|6^) with this asymptotic 
behavior. Power counting of this equation with respect to 
<fi gives then the roughness exponent £ RF = (4 — D)/3 
independent of N. But we are interested also in the fixed 
point function R*{4>) and the RG flow toward it. Using 
the ansatz R(<p) = R(\(f>\), the flow equation reduces after 
R^ R to 



dR 



1, 



= (e - 4C)«(0) + C<W) + ~[R"m 2 - R"(0)R"(<t>) 



+ (A-1)<™^''(0) 



R'{<j>) 



(64) 



The fixed point solution R*(4>) of this equation can be 
parameterized by its curvature A* = —R* (0), which is a 
measure for the disorder strength, and by its characteristic 
length scale £*. With these parameters, the fixed point 
function can be written as 



ir(0) = A*c 2 r(HC)- 



(65) 



This relation defines the dimensionless function r(u) with 
u = (/)/£,* and r"(0) = —1. It describes the functional form 
of the fixed point solution and has a characteristic scale 
of order unity. Notice that there is a hole family of fixed 
points which are parameterized by different £* since one 
has the freedom to choose an arbitrary overall length scale 
for <fi in the Hamiltonian ([j]) without crystal potential 
Vl (</>)• As we will see below, once a particular value for £* 
has been chosen, the corresponding A* is automatically 
fixed. From Eqs. (^34|) and ( |6^ ) follows that the function 
r(u) has to fulfill the equation 



= ar(u) + (3ur'(u) + ^[r"(u)} 2 + r" (u) 
Jl[r»] 2 r'(u) 



(66) 



where a = (e — 4()£* 2 /Z\* and [3 — (£,* 2 /A* are numerical 
coefficients. From this one obtains the exponent 



c 



(67) 



Generally, these coefficients are fixed first by the condition 
that r(u) has the correct asymptotic behavior and second 
by the choice of an overall length scale for <fi. From r(u) ~ 

— \u\ for large u follows the relation a = —f3. We choose 
a scale for <p such that r(0) = — 1. That r(0) has to be 
negative can be seen as follows. Evaluating Eq. ( |66] ) for 
small u, one easily obtains a — N/2r(0) due to r^O) = 

— 1. Since (3 = ££* 2 / A* > we have a < and our choice 
of scale corresponds to a = —N/2. 

After we have determined the coefficients a and /3, the 
fixed point function can be calculated in principle at least 
numerically by solving Eq. (|6^) for r(u) with initial condi- 
tions r(0) = — 1 and r'(0) = 0. But from a physical point 
of view, we are also interested in the evolution of the dis- 
order under the RG transformation sufficiently close to 
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the fixed point to study its stability or to calculate crit- 
ical exponents of an unstable fixed point. In particular, 
this becomes important if we take into account a compet- 
ing crystal potential which may lead to a new unstable 
fixed point associated with a roughening transition. Now, 
there is a special set of bare correlators Rq(4>) for which 
the RG flow can be calculated on all length scales exactly. 
This special set is given by the two parameter family of 
functions Ro(4») = A Qr(<f)/l;o), i.e., the functions are al- 
ready of the functional fixed point form but differ by their 
bare strength Aq and bare characteristic width £o- The 
RG flow in the hyper-plane which is formed by these two 
parameters can be obtained by inserting the Ansatz 



R(ct>) = A(l)e(l)rWai)) 



(68) 



into the flow equation (|64|). Using the fact that r(u) satis- 
fies Eq. d66| ) , this results in the following differential equa- 
tion for r (it), 



A(l)r(u) + B(l)ur'(u) = 0, 



(69) 



where 



A(l) = 2£A 
B(l) = -(A 



dl 
d£ 

dl 



,dA 



^■^-(e-AQeA + aA 2 , (70) 

(71) 



dl 

C^ 2 A + (3 A 2 



In the last formulas we have dropped the I dependence of £ 
and A. All the solutions r(u) of Eq. ( |69| ) are simple powers 
of u and vanish for u = 0. But this is in contradiction to 
the correct form of r(u) determined by Eq. (|6^), which is 
not a simple power of u and moreover finite for u = 0. 
Therefore such power law solutions cannot occur, and Eq. 
(||) is fulfilled if and only if A(l) = B(l) = for all I. This 
condition straightforwardly leads to the RG flow equations 



dA 
~~dT 

di 



(e-2QA 



< + f j 72 



a + 2(3 



A 2 



(72) 
(73) 



Flow equations of the same form but with different coef- 
ficients have been derived by Nattermann and Leschhorn 
[firjf for random bond disorder, using the approximation 
that the bare Gaussian shape of R(4>) is preserved under 
the RG flow. 

Before we come to the fixed point analysis of these 
equations, we will comment on the advantage of the two 
parameter representation over the functional flow of Eq. 
(p4j). For the Hamiltonian ([!]), without any perturbation 
from the crystal lattice, it has been shown that there exists 
a family of random field fixed points which is at least lo- 
cally stable and represents the rough phase of the interface 
p5| . ^3[ . Since this family of fixed points lies in the hyper- 
plane of the two parameter flow, the RG flow starting from 
an arbitrary bare random field disorder correlator Ro((f>) 
should be described on sufficiently large length scales also 
byEqs. ©, @. 



Let us now turn back to the fixed point analysis of Eqs. 
_T2]), ([73|). As explained above, we expect a line of stable 
fixed points resulting from the choice of the overall length 
scale for <p. This line is given by 



where 



A* 



C =£0 



a + 2(3 



4/3 A 



CA 



(74) 



(75) 



depends on the bare parameters £0 &ud Aq. The corre- 
sponding RG flow is shown in Fig. ^. Linearizing about 
the fixed points, one finds the universal eigenvalues Ai = 
and A2 = — e corresponding to a redundant marginal and a 
stable direction, respectively. The leading irrelevant eigen- 
value of the full functional RG flow for random field disor- 
der is indeed given by A2 [Q. But notice that in general 
one cannot expect to obtain by the two parameter ap- 
proach the leading irrelevant eigenvalue out of the infinite 
number of irrelevant eigenvalues. 




Fig. 5. Schematic RG flow to a line of fixed points in the 
hyper-plane formed by Aq and £o- 



The transversal displacement correlation function 
1 



C(x) = -<[0(x) 



0(O)F 



(76) 



can be calculated perturbatively from the renormalized 
but non-rescaled effective parameter A{1) = A{l)e ty2i '~'^ 1 
using I = ln(yl/g). On length scales larger than the Larkin 
length, |x| > L A = [eN/(a + 4/3)^/Z\ ] 1/e , we obtain 



C(x) = {N/(3)eM/L A ) 



2< 



(77) 



where we have explicitly considered the iV-depence of £0 
by writing £q = A£q. Notice that the Larkin length La 
and the coefficient of C (x) remain finite for N — > 00 due 
to a = -P = -N/2. 

For D = 4, the roughness exponent ( vanishes and 
logarithmic scaling of the transversal displacement is ex- 
pected. From the RG flow of the effective parameter A 
for e = 0, we find that the roughness grows only sub- 
logarithmically, 



CM 



ln° 



l x l/ a oJ 



2(3 



4(3 



(78) 
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For random field disorder, we have the exponent a = 2/3 
independent of N. 



5.1.2 Random Bonds 



As for random field disorder, it is straightforward to see 
that the functional behavior of the exponentially decaying 
correlator Rq(4>) for random bond disorder is preserved by 
the RG flow p3[ . To obtain a valid fixed-point function 
for short-range correlated disorder, therefore, we have to 
search for solutions of Eq. (|64|) with an exponential tail at 
large <f>. The fixed-point function can still be written in the 
form of Eq. (65), and it remains the problem to determine 
the coefficients a and (3 such that Eq. (^6|) has an expo- 
nentially decaying solution r(u) for the initial conditions 
r(0) = 1, r'(0) = 0. As in the random field case we have 
a = N/2r(0) = N/2 by our choice of scale. Unfortunately, 
the conditions for the behavior of r(u) do not fix the co- 
efficient (3 uniquely since there exists a discrete family of 
solutions of Eq. ( pq ) which decay exponentially at infinity. 
But it can be shown that only the function of the whole 
set of solutions corresponding to the largest value of (3 is 
stable with respect to short-range perturbations Q . The 
corresponding value for [3 has to be calculated numeri- 
cally. Solving Eq. (^) numerically, the unstable solutions 
can be distinguished from the stable one since they show 
oscillations instead of a simple exponential decay. The re- 
sulting values of (3 are shown in Tab. ^ together with the 
corresponding roughness exponents C for different N. The 
values for (3 decrease with N and approach 1 /2 for N — > oo 
corresponding to £/e — > 1/(4 + N). 



N 


P 


C/e= (4 + a//3)- 1 


1 


0.6244 


0.2083 


2 


0.6010 


0.1766 


3 


0.5807 


0.1519 


4 


0.5635 


0.1325 


5 


0.5491 


0.1169 



Table 2. Numerical values for /3 and the exponent 



With these values for the coefficients a and (3 the flow 
equations Eqs. (f72|) and (^3|) remain valid also for the 
random bond case. The RG flow is still given by Eq. ( [74] ) 
but with the exponent C, from Tab. 2[ and looks similar to 
that shown in Fig. ||. For N < oo, the correlation function 
C(x) is given by Eq. (|77]) but with the values for £ from 
Tab. |[ In contrast, for N — oo the roughness exponent 
vanishes and the manifold is only logarithmically rough, 



C(x)=4 e | 2 ln(|x|/i zl ). 



(79) 



The exponent a describing the sub-logarithmic scaling 
of Eq. © valid for D = 4 becomes a = 0.4166 for N = 1 
and approaches a = 2/(N + 4) for N — > oo. 



5.2 Periodic Elastic Media 

The periodicity of the unrenormalized Rq(4>) has to be 
preserved under the RG flow, hence we have to search 
for fixed point solutions of Eq. (|6^) with an appropriate 
symmetry given by the reciprocal lattice vectors Q. In 
general, one has to treat different lattice symmetries sep- 
arately. For the simplest case, a square lattice, the fixed 
point solution is separable, 



TV 



(80) 



where R*(<f>) is a periodic fixed point function of Eq. (J6S 
with N = 1 and the 0™ are the components of <fi. For 
a triangular lattice the situation is more difficult and we 
did not try to find a solution for this case. The N = 1 
fixed point solution R(cj>) corresponds to the charge den- 
sity wave case and has been calculated in Ref. [B9|. The 
RG flow towards this fixed point can be analyzed by keep- 
ing track of the flow of the Fourier coefficients of R(<fr). 
For simplicity, we choose the periodicity as 2tt and take 
R{4>) — > R{4>)- Inserting the Fourier expansion 



R((j)) = ^ R m cos(m0) 



(81) 



■m— 1 



into the functional flow Eq. ( |62] ) where we have to choose 
C = due to the periodicity, we obtain the following in- 
finitely many flow equations for the coefficients R m which 
are coupled to each other, 



dR n 



771 — 1 

eR m + ^ X! m ' 2 ( m ~ rn') 2 R m ,R n 

rn' — l 

^ oo 

- ^ m ' 2 i m ' - m ) 2 R m' Rm' -m 



m'— m+1 
„ , 2 



(82) 



By performing simply the sums, it can be checked that 
a fixed point of this set of equations is given by R^ = 
2e/3TO 4 . The corresponding series (81) can be summed up 
and yields the solution known from Ref. |39|. The flow 
towards this fixed point as well as its stability can be an- 
alyzed by rewriting the Fourier coefficients as 



Rr, 



6 A 



+ r v 



(83) 



where A = —R"(0) with the fixed point value A* = ctt 2 /9. 
Substituting this Ansatz into Eq. ( |S2| ) , we see that the r m 
flow to zero with eigenvalues A m = e(2/3 — m — 2m 2 /3) 
whereas the flow of A is determined by 



dA 
If 



eA 



9 



A'. 



(84) 
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Since we are interested in the RG flow on sufficiently large 
length scales, we will therefore assume here the particu- 
larly simple functional form of the fixed point by dropping 
the r m . This approximation is justified on length scales 
larger than the Fukuyama-Lee length La = A^ 1 ^ where 
the parameters r m approach zero. The RG flow towards 
the fixed point is then described by the single parameter 
flow of Eq. 



After projection onto the lowest harmonic, this equation 
reduces to 

1 . ,(2tt) 6 



V [W((f>)] 



-An* 



since sin 2 (y) = (1 — cos(2j/))/2. Note that the negative 
eigenvalue 9 of the temperature flow, see Eq. (^lj), is also 
shifted by this constant to a larger value. Physically, this 
corresponds to an increased elastic stiffness of the mani- 



The transversal displacement correlation function growthsf ld due to the lattice potential. 



only logarithmically on scales beyond the Fukuyama-Lee 
length La, 



C(x) 



2— eln(|x|/i 4 ) 



(85) 



with a universal coefficient which may depend on the ac- 
tual lattice symmetry. The coefficient given here applies 
to a square lattice with a periodicity of 2n. 



6 The roughening transition 

To study the effect of a finite lattice pinning potential 
Vl{4>) upon the large scale behavior of elastic objects, we 
have to analyze the full set of functional RG equations 
(|6l|)-(|63|). As has been shown is Section 5, while the same 
functional RG equations apply for the different types of 
elastic objects and bare disorder correlations, the corre- 
sponding solutions are quite different. Therefore, we will 
discuss the analysis of the RG flow in the following subsec- 
tions separately for elastic manifolds and elastic periodic 
media. 



6.1 Elastic manifolds 

For the case of elastic manifolds, the functional RG flow 
of the disorder correlator R(4>) has been reduced already 
to a flow of two parameters for vanishing lattice pinning 
in Section 5. To develop the RG analysis for finite Vl(<P), 
it is important to note that the functional form of the flow 
equation for R(4>) is not changed. Indeed, a finite lattice 
potential leads only to a shift of e = 4 — D by the constant 
2Vq[W(4>)], as can be seen easily from Eq. (p2[). Therefore, 
the RG flow of R(cf>) can be replaced again by the two- 
parameter flow of Eqs. ( [72|) and ( |73|) but with a shifted 
e. For the lattice potential we have to consider only the 
first harmonics, and we assume in the following the simple 
form of Eq. (0) with 



V{4>) = v cos ( — (j) 



(86) 



for all components of the vector <fi. 

To begin with, we calculate the projection of the func- 
tion W{(j>) defined in E q. (]60| ) onto the lowest harmonic 
E {4>) = 1. Using Eq. @ and d l d J R(0) = -ASy, the 
function W{<p) becomes 



W{<j>) 



-Av 



2 (2tt) 



6 N /2tt 



(87) 



The RG flow of the first harmonic of the lattice poten- 
tial can be obtained by insert ing the sum of the first har- 
monics ( j86] ) for Vl(4>) in Eq. ( |63| ) and projecting again the 
whole equation onto the first harmonics. To treat the sec- 
ond term ~ (j? in Eq. ([53J) properly, one has to rescale also 
the periodicity a which leads to the trivial flow da/dl = 
—C,a. After the projection, the term which is quadratic in 
Vl vanishes since it contributes only to the second har- 
monics. Both the term cubic in Vl and the term 7 7 o[W] Vl 
from the stiffness renormalization are of order v 3 . Using 
for the projection of W the relation cos 3 (j/) = (3cos(y) + 
cos(3y))/4, we obtain the following parameter RG flow, 



-7T = (^-2C)^ 



dT 
dl 
dA 
~Jl 

d A 
di 

dv 

~dl 

da 
~dl 



2-D-2 C - 1 -Av 2 ^ 
2 a 6 



T. 



a + 2/3 (2tt) 6 



2 - 2C - - 
2 

-Co. 



1 (27T) 2 



A 



7 (2ir) e 
v 



8 a b 



(89) 



(90) 



(91) 



(92) 
(93) 



It is important to note that these flow equations are sym- 
metric with respect to v — > —v as could be expected for 
a periodic potential since this transformation corresponds 
only to a shift in <fi l by a/2. Naively, one expects for the 
coefficients a and (3 the same relations as given below Eq. 
( |66| ) but with shifted e. But these relations have to be re- 
examined later since it turns out that £* 2 /A* diverges at 
the new fixed point describing the roughening transition. 

To analyze the above scaling transformations, it is use- 
ful to introduce the new variables A = (2n/a) 2 A, £ = 
(27r/a)£ and v = (2n/a) 2 v. By this change of variables, 
the roughness exponent ( is eliminated from the coupled 
flow equations. To finally obtain £, we have to evaluate the 
relations for a or (3 in terms of the flow of the new vari- 
ables along the critical line flowing into the fixed point. In 
terms of the new variables, we have 



dT 
~dl 
dA 
~dT 
d? 
dl 
dv 
~dl 



2-D-2C- -Av 2 
2 



eA- 

2(3A, 



a + 2(3 



T. 

A 2 , 



2- -A 
2 



7 A 3 
8 



(94) 

(95) 

(96) 
(97) 
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Before we study the case v 7^ 0, it should be mentioned 
that the stable fixed point of the rough phase is located 
at v — 0, A — 00 in terms of the new variables since the 
periodicity a flows exponentially fast to zero. Now we turn 
to the fixed point analysis for finite v and N < 00. The 
case N = 00 will be discussed separately in Section 6.4. 
To obtain a fixed point with a finite value for v from Eq. 
(|97|) , the parameter A has to flow along a critical line to 
a finite value A* . The corresponding fixed point value for 
v is then given by v* 2 — 4(4/ A* — l)/7. As a consequence 
of Eq. (96) and since (3 > 0, one obtains the asymptotic 
behavior^ 2 = 2(3 A* I up to sub-linear corrections. But this 
divergence of £ 2 means that we can drop the first term in 
the brackets of Eq. ( |95|) for the fixed point analysis. This 
leads to the fixed point values, valid to lowest order in e, 



A*=A-\e, 



c 

T 



(98) 



As one may expect, v* vanishes when the upper criti- 
cal dimension is approached, justifying the e-expansion. 
It should be noted that v* has the meaning of a (rescaled) 
mass for the field </>. The second small parameter in the 
double expansion with respect to both pinning potentials 
is A, which is indeed small at the new fixed point since 
A ~ l/l along the critical line as we will sec below. More- 
over, it should be mentioned that the fixed point value A* 
vanishes for e = 16/7 corresponding to a shift of this fixed 
point to the v axis for D = 12/7, see Fig. |^. Note that 
this is, despite the lowest order e result, in approximate 
agreement with the absence of a roughening transition for 
D < 2 as predicted by scaling arguments in Section 3. 

Linearizing around this fixed point, we obtain, in low- 
est order in e, the universal eigenvalues A± = ±2^/e. Since 
the fixed point has an unstable direction, it has to be as- 
sociated with the roughening transition. In the flat phase, 
the size of typical excursions of the manifold from the pre- 
ferred minimum of the lattice potential defines the longi- 
tudinal correlation length 



1 



(99) 



Interestingly, m does not dependent on a, (3 and is there- 
fore universal for random field and random bond disorder. 
The same expression for v\\ has been obtained previously 
for the thermal roughening transition but with e = 2 — D 

§3- 

To determine the critical roughness exponent £ c at the 
roughening transition, we have to re-examine the relations 
between the coefficients a and (3 and the fixed point values. 
Paying attention to the fact that A approaches zero at the 
new fixed point, we obtain the refined expression 



r ? 

a = urn -=■ 

l^oo A 



a 2 1 dA 

2 ^ A "-A^T 



which turns out to be equivalent to 
(3= lim CV 



(100) 



(101) 



From the last expression we see that £ has to flow along 
the critical line asymptotically according to 



c 



1 

21 



(102) 



since £ 2 / 'A = 2(31. Therefore, at the transition we have 
( c = and the roughness of the manifold grows only log- 
arithmically. Due to ( = 1/21 for large the periodicity 
of the lattice potential flows only algebraically to zero, 
a ~ l^ 1 / 2 , in agreement with the logarithmic scaling at 
the transition. Therefore, the renormalized width £ ~ £/a 
of the disorder correlator flows to a finite value £* at the 
new fixed point. Notice that at this fixed point one has no 
longer the freedom to choose an arbitrary overall length 
scale for (p due to the fixed periodicity ao of the lattice 
potential. As a consequence, there is no line of fixed points 
connected by a redundant operator as in the case of van- 
ishing lattice potential. 

The free energy fluctuations grow with length as L e , 
where 8 is defined by Eq. (|6l|). In the rough phase the 
tilt symmetry of the Hamiltonian (l|) guarantees that the 
temperature is renormalized only by scale changes. This 
leads to the exact result 9 = D — 2 + 2(. The relevance of 
the lattice potential results in an additional renormaliza- 
tion of temperature and therefore a different exponent 8 C 
at the transition. Independent of a and (3, we obtain from 
the fixed point values of Eq. (^8|) to lowest order in e the 
result 

0c = 2--. (103) 

The cross-over length scale Lr introduced in Section 3 
can also be obtained from the RG flow of A. This length 
corresponds to the scale where the roughness of the mani- 
fold changes from algebraic to logarithmic growth on both 
sides of the transition. In terms of the RG flow, Lr = 
aoe lR is determined by the scale Ir_ where A(l) approaches 
its asymptotic fixed point value A*. To estimate Ir, one 
can integrate Eq. (95) with v = to obtain a primary 
solution. The effect of the w 2 -term can subsequently be 
included by matching the primary solution A{1) and A* 
at Ir. Doing so, we get 



U 



■ 4/3 A* a 2 



1/2C 



N 



(2 
Ml 



(104) 



where £ is the roughness exponent of the algebraic regime 
at smaller scales. It is given by the exponent £ of the 
rough phase for £ "C ao an d C = e /2 for £ ^ ao- This 
is in agreement with the scaling result of Eq. (|l^) if vq is 
replaced there by w ,c of Eq. (|l7|). 

Beyond the length Lr and sufficiently close to the 
transition where £11 > |x|, the asymptotic behavior of the 
transverse difference correlation function is given by 



C(x) = -<[0(x) - 0(O)H ~ ln(|x|/L«) (105) 

on both sides of the RT, see Fig. ||(d). Beyond £11 the 
roughness crosses over to the power law C(x) ~ (|x|/£||) 2 '' 
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in the rough phase. On the flat side of the transition 
C(x) saturates on scales larger than £|| at a finite value 

The phase boundary between the rough and the flat 
phase can be obtained from the condition that the primary 
solution of Eq. (|97]) , neglecting the v 3 term, and the fixed 
point value v* match at the cross-over length Ir. This 
condition leads to the result 



2/C 



i>o, c (A},£o) = g(e)a — 3- 



a 




2/e 



(106) 



which is in agreement with the estimate of Eq. (|2J) de- 
rived from scaling arguments for random field disorder. 
The coefficient g(e) is given to lowest order in e by 



9(e) 



eN 



4/3 



l/C-2/e 



(4e) 



(107) 



In both expressions above, C is always the exponent of the 
rough phase for £ <C ao and Q — e/2 for £ ^> a o- As 
can be seen easily from the behavior of 17(e), the thresh- 
old i>o, c vanishes exponentially fast for e — > 0. This is in 
agreement with the fact that in four dimensions an arbi- 
trarily small lattice potential leads to a flat phase due to 
the sub- logarithmic roughness for v = 0, cf. Eq. ([78]). The 
schematic RG flow for 2 < D < 4 is shown together with 
that expected from the scaling arguments of Section 3 for 
D < 2 in Fig. §. 



6.2 Interfaces in Dipolar Systems 

As shown in Section 3, in the physically interesting case 
D = 2 the manifold is rough beyond the cross-over length 
Lr even for a finite crystal potential. Although Lr is ex- 
ponentially large for weak disorder, there exists, strictly 
speaking, no roughening transition for D — 2. However, 
the transition is expected to be seen even for D = 2 in 
systems with dipolar interactions as already noted in Ref. 
pq|. In these systems ferroelectric or magnetic domain 
walls (N — 1) separate regions of opposite polarization or 
magnetization. Assuming that the dipoles interact via a 
three-dimensional Coulomb force and that the dipolar axis 
p is parallel to the wall, one obtains for D > 1 in harmonic 
approximation the modified elastic energy |47[f48|] 



E c \ = 




(108) 



where g measures the relative strength of the dipolar inter- 
action. As can be easily seen from Eq. ( |108| ), the dipolar 
interaction increases the elastic stiffness of the interface 
on large scales. As a result, the upper critical dimension 
D c should decrease compared to the case of an isotropic 
elastic kernel ~ q 2 with D c = 4. 

Indeed, our results can also b e ap plied to interfaces 
with the elastic interaction of Eq. ( 108 ). For the RG anal- 
ysis of dipolar interfaces, there are mainly two differences 




Fig. 6. Schematic RG flow for TV < 00 as a function of the 
disorder variable A and the lattice strength v for (a) D < 2, (b) 
2 < D < 4 and (c) D > 4. Note that the fixed point associated 
with the rough phase is shifted here to A — 00. (d) shows 
the length scale L dependent roughness near the roughening 
transition where the correlation length £y (v) diverges. 



compared to the usual case discussed before. First, the 
large scale behavior of the interface is described correctly if 
we take the stiffness constant to be an effective q-dependent 
function 7(q) ~ y/gj]^ 1 ^ ^ 2 ■ This can be seen by writ- 
ing q p — q cos(iy9) and performing the angular integra- 
tion over ip in the Fourier representation of the free two- 
point function. This integration yields an additional factor 
~ Iql' 15-1 )/ 2 in the limit |q| — > which can be absorbed 
in 7. Repeating the RG analysis with this modified 7, the 
form of the terms in the flow equations arising from scale 
changes remains unchanged if we replace e = 4 — D by 



(109) 



The second difference arises from the momentum shell in- 
tegral over the anisotropic elastic kernel. This leads only 
to a modified numerical coefficient Kd in the non-linear 
parts of the flow equations, but has no influence on the 
exponents. Therefore, the exponents calculated abo ve re- 
main valid to first order in e if we use Eq. ( 109 ). The 
upper critical dimension is shifted to D c = 3, and there 
exists a roughening transition for two dimensional dipolar 
interfaces which is described by our results with e = 3/2. 



6.3 Periodic Elastic Media 

Having established the existence of a roughening transi- 
tion for single elastic manifolds, we now turn to the case 
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of periodic elastic media. From the analysis of Section 5.2 
we know that the correlation function C(x) behaves log- 
arithmically for |x| > L/\ in the rough unlocked phase 
(U) described by the stable T = fixed point U*. There- 
fore, one can expect a logarithmic scaling of <f> also at the 
roughening transition, hence we set in this section £ = 
from the very beginning. To extend the RG analysis to 
Vl{4>) 0; we follow the treatment presented above for 
elastic manifolds. In the following, we will focus on the 
particular case of a square lattice for the elastic medium. 
Then we are able to make use of the reduction of the 
i?(0)-flow to the single parameter flow for A derived in 
Section 5.2, cf. Eq. (pif). Again, a finite lattice potential 
changes this equation only by a shift in e. For simplicity, 
we set the periodicity of the elastic medium to I — 2n 
in the following such that p = 2n/ao on all length scales 
since transversal lengths are not rescaled here. Then the 
scaling equations for elastic media read 

^ = {2-D-\p G Av 2 )T, (110) 

^ = eA-^A 2 -p 6 v 2 A 2 , (111) 
al tt z 

From the scaling arguments presented in Section 3, one 
may expect that for p < p c — weak disorder is 

always irrelevant due to the logarithmic roughness in the 
phase U. Forp > p c a weak periodic potential is irrelevant, 
but one expects that it becomes a relevant perturbation 
once vq exceeds a critical value vq_ c . Therefore, vo, c h as to 
vanish if p — > pj, and a new unstable fixed point corre- 
sponding to a roughening transition should be accessible 
perturbatively in e = A—d and the second small parameter 
H = p 2 /p 2 c -l. 

Indeed, there exists a new unstable T — fixed point 
F* which approaches the stable fixed point describing the 
rough phase if /i — > 0. It is parameterized to lowest order 
in e and \i by 

A* = ^e(l-»), p 2 v* = ^. (113) 

Linearization around this fixed point gives the eigenval- 
ues Ai = — e and A2 = 4/i. Thus, one has to identify this 
unstable fixed point with the roughening transition. The 
longitudinal correlation length £|| beyond which the pe- 
riodic medium is locked into the minima of the lattice 
potential diverges according to 

Z\\~\v-v , c \-"", i/|| = i (114) 

if the transition is approached from the flat side. On the 
rough side of the transition £|| is infinite, in contrast to the 
situation for single elastic manifolds. The resulting phase 
diagram is characterized by a critical line A c (v) separating 
the flat or commensurate (C) and rough, unlocked phase 
(U), see Fig. 0. The position of this critical line is given 
by the scaling result i>o, c ~ P 2 ^o derived in Section 3. 




1 1+ V 



Fig. 7. Schematic RG flow for p > p c in the w-Zi-plane. A 
critical line A c (v) with fixed point F* separates the flat (G) 
and rough, unlocked ( U) phase. For p < p c the fixed points F* 
and U* merge. Since for large disorder strength A topological 
defects will proliferate, we expect a disordered phase D above 
the line Ad(v) as will be explained in Section 8. While the RG 
flow in this range is not yet clear, the transition from the G to 
the D phase is probably in the universality class of the random 
field p-state clock model. 

To characterize the order of the periodic medium, we 
introduce the translational order parameter 

!P(x)=piexp(*0(x)). (115) 

For simplicity, we study here the case of a scalar field (N = 
1) only. Inside the C phase, one has a nonzero expectation 
value (^(x)). The correlation function of the fluctuations 
5&(x.) = ^(x) — (^(x)) decays exponentially in this phase, 

if(x) = (^(x)^*(0)) - exp(-|x|/£||). (116) 

Approaching the transition from the C phase, the order 
parameter vanishes continuously as 

Wx»~(vo-V 0tC )P, /?=^1/||. (117) 

In contrast, at the transition and in the rough phase, vq < 
i»o,o the periodic medium is unlocked with (^(x)) = and 
the correlation function shows a power law decay, 

K(x) ~ \^- D -\ fj= (l + 7r 2 /9)e. (118) 

At the transition the 'violation of hyperscaling' exponent 
is increased from the exact value 9 = D — 2 in the rough 
phase to 

6 C = D -2 + efi/2 (119) 

due to the renormalization of the elastic stiffness of the 
medium. The critical exponents mentioned above fulfill 
all scaling relations and exponent inequalities which hold 
for random field models with a discrete symmetry of the 
order parameter ]49| , |50|| . 
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6.4 The Limit N = oo 

It is important to note that the order as well as the critical 
exponents uu and £ c of the roughening transition do not 
depend on TV as long as N is finite. For arbitrary N, the 
roughening transition of manifolds in random bond sys- 
tems has also been studied within a variational approach 
by Bouchaud and Georges (BG) ||. At T = and for 
£o ao these authors find a first order transition between 
a glassy rough and a glassy flat phase. In the glassy rough 
phase they obtain the Flory result £ = e/5 which is known 
to be not correct for small e. We attribute this as well as 
the first order nature of the transition to the use of the 
variational method which is known to give spurious first 
order transitions |5l|| . However, the variational approach 
is expected to become exact for N — oo. Since our results 
seem to be in contrast to the N — oo results of BG, we 
will study the RG flow for elastic manifolds in this limit 
now. 

Our discussion starts again with the flow Eqs. (|95|)- 
( p7| ) . To have a finite bare starting value for £ in the limit 
N — > oo, we rescale this parameter according to £ 2 — > 
N£ 2 . For random field disorder, the structure of the flow 
equations remains unchanged in the limit N — > oo since 
(3 = N/2. Therefore, the transition is still o f sec ond order 
and the phase boundary is given by E q. ( 106 ) with the 
N- independent coefficient of Eq. (107). In contrast, for 
random bond disorder the manifold is only logarithmically 
rough in the case N = oo, cf. Eq. (|79|), and we expect a 
roughening transition similar to that of periodic media. 
Indeed, due to the above rescaling of £ and — 1/2 we 
obtain d£ 2 jdl — and, therefore, Eq. (|95|) reduces to 



dA 

~dT 



i 

V 



A 2 



(120) 



where p = 2-7r£o/ao is a con stant. Comparing this equa- 
tion with the flow Eq. (Ill) of A for periodic media, we 
conclude that for p < p c = y/2/e the manifold is already 
in the flat phase for an infinitcsimally small vq. A rough- 
ening transition at a finite vq exists only for p > p Cl and 
t hen the correlation length diverges as described by Eq. 



114) 



The variational calculation of BG applies in the ran- 
dom bond case where £o -C ao, i.e. for p < p c . Therefore, 
for N — oo we expect no transition at a finite Vq. Indeed, 
this is in agreement with the result for the phase bound- 
ary of BG in the limit N = oo. The resulting transition 
at vo, c = for p < p c is not described by a fixed point of 
our RG equations and we have been unable to ascertain 
definitely the nature of the transition in this case. We will 
discuss this in more detail in Section 8 in the context of 
periodic media. 



6.5 Incommensurate Phases 

So far we have considered only periodic elastic media with 
a periodicity I = pao which is exactly an integer multiple 



of the crystal lattice periodicity ao- But in general, this 
condition is not fulfilled exactly in real systems and one 
has to assume a shifted lattice constant 



I =pa (l -pa S) 



(121) 



for the elastic medium, where S is a small misfit. For ex- 
ample, 5 can be changed in flux line lattices by tuning 
the external magnetic field or in charge density waves by 
changing the temperature. It is well known from the sine- 
Gordon model subject to thermal fluctuations that above 
some critical value S c of the misfit, the incommensurabil- 
ity is compensated by the formation of domain walls or 
solitons of intrinsic width given by the correlation length 
£|l ]52fl . Across these walls the phase field 0(x) changes 
by 2n/p. In the following, we specialize our discussion of 
the incommensurate phases to the case of a scalar field 
(N = l). 

The effect of a finite misfit can be described by adding 
a linear gradient term to the Hamiltonian defined in Eq. 
(pi), such that the elastic term is replaced by 



J d D x{V(f>-5} 2 



(122) 



with S — Sz and z the modulation direction of the periodic 
medium. For completeness we note that for a vanishing 
lattice potential Vl the misfit 5 can be eliminated by the 
transformation </>(x) — > <fi(x) + 5z due to the statistical 
symmetry of the diso rder under a shift in <f>. The linear 
gradient term in Eq. (122) can be integrated out leading 
to the contribution 



(123) 



to the Hamiltonian (0), where we have introduced the 
mean soliton density i -1 — p[<fi(z = L) — <fi(z = 0)]/27rL 
for a system of size L in z-direction. 

To study the case of finite lattice pinning, a direct 
RG treatment of the Hamiltonian (0) with an additional 



term (123) is not justified since a finite soliton density l~ x 
induces linearly growing contribution ~ £~ 1 z to the field 
</>(x) which cannot be Fourier expanded. To avoid this 
problem, we go back to periodic boundary conditions by 
the transformation </>(x) — > </>(x) — 2it z/pl |5j|. Note that 
we treat the mean soliton distance i as a fixed parameter 
which is not renormalized and has to be determined later 
by minimizing the free energy of the soliton lattice. In the 
commensurate phase (C) for 5 < S c , the soliton density 
vanishes and hence the RG analysis for 6 — applies. 
To estimate 5 C , we compare the energy gain of a soliton 
per unit area, — 7<527r/p, with its cost for vanishing misfit, 
Jy/Vcff^/p, where v c s — (p£||) -2 denotes the effective 
strength of the lattice potential renormalized by disorder. 
This leads to the result 



Pt\\ 



VQ 1 



An 



J 0,c 



(124) 



where we have used Eq. (114). The corresponding phase 
diagram for fixed Vq is shown in Fig. |[ Since for p < p c 
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the lattice potential is always relevant for weak disorder, 
there exists in this case for large enough disorder a di- 
rect transition to a disordered phase (D). This transition 
is not described by our approach since it is driven by the 
proliferation of topological defect s as will be explained in 
Section 8. In deriving expression (124), we have neglected 
a possible additional reduction of 8 C by fluctuations of the 
solitons on scales beyond £|| which have not been inte- 
grated out by the RG. 




(b) 




Fig. 8. Schematic phase diagrams for finite misfit 5 at fixed 
crystal potential strength vo for p > p c (a) and p < p c (b). Due 
to the proliferation of topological defects a disordered phase is 
expected above a disorder strength Ad ■ 



In the incommensurate phase (I) the above transfor- 
mation to periodic boundary conditions induces a spatial 
oscillation in the amplitude of the lattice potential, 



Vl(</>) = vq cos(p0 — 2irz/£). 



(125) 



Notice that the same oscillating pinning potential describes 
vicinal surfaces which are tilted by a small angle 8 — 2tt/p£ 
with respect to the crystal planes Q. In the RG calcu- 
lations, this leads to an additional factor exp(i2irz/£) in 
the integral of Eq. ([l|]) . When the lattice potential is rel- 
evant, Aq < Aq >c , there are two length scales: £ and the 
correlation length £|| . 

If £|| <§C £, the RG up to £|| is insensitive to the oscilla- 
tions along the z-dircction and the soliton lattice can be 
described in terms of an effective free energy. The corre- 



sponding free energy density is given by |48 



all 



4cre 



B(D)ai 



-IIP 



(126) 

Here the first term is the free energy of a soliton per 
unit area, the second term is the bare soliton repulsion, 
and the last part describes the steric repulsion caused by 
collisions of the solitons. The effective soliton tension is 
a = 8j/p 2 £\\, B(D) is a coefficient which vanishes for 
D > 5 and /? = </(2 - 2Q = (5 - D)/2(D - 2) where 
£ = (5 — D)/3 is the roughness exponent of the soliton 
walls subject to random field disorder. Minimizing the free 
energy, we obtain the mean soliton distance near the IC- 
transition 

e~(6-5 B )-P. (127) 



To determine the translational order inside the I phase, 
the soliton lattice itself can be considered, on scales larger 
than £ as a periodic elastic medium with a disorder corre- 
lator Ro(4>) of periodicity I, The corresponding Hamilto- 
nian is given by Eq. ([!]) with Vl = and N — 1 but with 
anisotropic elasticity. The elastic constants are 7|| = a/£ 
parallel to the solitons and 7j_ = £ 2 f"(£) in z direction. 
Rescaling the z coordinate according to •v/7||/7±z — > z, 
we can achieve again an isotropic elastic term with 7 = 
^711 7j_. Using the result of Eq. ( |85| ) we obtain a loga- 
rithmically diverging correlation function for the soliton 
displacement field it(x), 



([u(x)-u(0)P) = 



A-D 

i 

36 




(128) 



Here the Larkin lengths L\\ ~ (i 4 ^ 2 ^ 2 ) 1 /^-^ and 

■k-L = vTjlTtP^N are the correlations length beyond which 
the logarithmic divergence becomes asymptotically exact. 
Since 7|| ~ t~ l and 1± ~ £K" 2 )/« = £-(n+i)/(5-D) for 
random fields, both correlation lengths diverge if the IC- 
transition is approached from the I phase as 

L|| „ ^/(5-D) ^ {s _ §c y 3 /2(D-2) (12g) 

L ± ~ £ ~ (5 - <y-C5-D)/2(D-2) (13Q) 

similar to the situation for the IC-transition in the pres- 
ence of thermal fluctuations only (4^^] . 

The behavior of the correlation function if(x) of the 
fluctuations of the translati onal order parameter ^(x) can 
be easily obtained from Eq. (128). The soliton lattice shows 
quasi long range order with 



^(x)=P?cos— 



-(4-£>)7i- 2 /18p 2 



(131) 



Here we have used that in the I phase a difference in the 
soliton displacement is related to a shift of the phase ^(x) 
of the order parameter by 0(x) — 0(0) = (2ir/p)(z — u(x) + 
u(0))/£. 

In the opposite case where ^|| 3> £, the lattice po- 
tential Vl should be irrelevant since it averages to zero 
on scales larger than £ for geometric reasons. This corre- 
sponds to the situation above the roughening transition, 
Aq > Aq )C , and now the original periodic medium (instead 
of the soliton lattice) is described by the Hamiltonian (0) 
with Vl = since the misfit can now be shifted away as 
mentioned above. Therefore, the correlation function has 
t o sh ow a crossover to a behavior which is given by Eq. 
(131) where the Larkin lengths of the soliton lattice be- 



come that of the original periodic medium, Ln , L ±_ — > 



The exponent of the algebraic decay in Eq. (131) goes over 
to -(4- L>)tt 2 /18. 



7 Experimental Implications 

Having established the possibility of a disorder driven 
roughening transition for various systems which can be 
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described in terms of elastic objects, we discuss now the 
experimental consequences of our results. To begin with, 
we comment on how the transition may be approached 
for a sample with fixed disorder strength. In general, for 
both elastic manifolds and periodic elastic media one may 
expect that the bare pinning strengths v and A de- 
pend also on temperature. Hence the transition should 
be reachable by changing the temperature T for appropri- 
ate T-dependencies of both pinning strengths. The pre- 
cise T-dependence of the parameters is model dependent 
and, therefore, we will give only some general arguments 
and typical examples. Close to the condensation temper- 
ature T c of the system where the intrinsic width of the 
elastic manifold becomes much larger than the periodicity 
of the lattice, £o 3> ao, we expect a very weak influence 
of the crystal potential. An example are domain walls in 
magnets where T c is the Curie temperature. For Peierls 
barriers which appear, e.g., for domain walls in ferroelas- 
tics or solitons in incommensurate phases, one expects an 
exponentially decreasing influence of the lattice potential, 
u ~ e~ c ^°. 

The most prominent example for periodic elastic me- 
dia are charge density waves. In this case the critical ra- 
tio of the pinning strengths depend also on temperature 
since vq/Aq ~ p\ +p for D = 3 where the charge density 
amplitude p\ changes with T and vanishes at the Peierls 
temperature. Thus, for not too large disorder, we can ex- 
pect to see the roughening transition by increasing T for 
both elastic manifolds and periodic media. 

Indications about which phase is present should be 
given by the creep behavior since at small driving forces 
the dynamics are governed by the static behavior of the 
elastic objects. Indeed, under the influence of a small driv- 
ing force density / ex , the motion of the elastic object 
is dominated by jumps between neighboring metastable 
states in the rough phase and between adjacent minima 
of the lattice potential in the flat phase. As has been dis- 
cussed earlier |2^Jil]j5^| , in both cases the creep velocity 
u is exponentially small, 



bond or random field disorder, 



w(/ex) ~ exp 



Ec 
T 



(132) 



where f c denotes the maximal pinning force density of the 
dominating pinning mechanism and E c the correspond- 
ing energy barrier. The creep exponent k depends on the 
phase which is present and on the type of disorder under 
consideration, i.e. on the disorder correlator Rq (<fi). 

For periodic elastic media, we obtain in the rough un- 
locked and flat phase, respectively, 



D-2 



trough 



Kflat 



D-l. 



(133) 



The pinning force density and the energy barrier in the 
rough phase are given by f c s» "fL^ and E c w jL®~ 2 , 
whereas in the commensurate flat phase they are f c ~ 
7 ^ 2 and£; c «7^- 2 . 

For elastic manifolds, the exponent n depends on the 
type of disorder via the roughness exponent Q , i.e. random 



trough 



D - 2 + 2C 
2-C ' 



Kflat 



D-l. 



(134) 



lHL°- 2 and f c 



In the rough phase, we obtain E c 

7£ £^ 2 , and in the flat phase E c ~ j V ( 2 - D )/ 2 and f c ~ jv 
for an interface with short range interactions. For dipolar 
systems one has to replace the dimension D according 
to Eq. (109), which gives k rf = 1, k rb = 0.6666 for a 
two dimensional interface in the rough phase and k = 1 
for a flat interface independent of the type of disorder. 
Therefore, in the RB case the phase of the interface can 
be determined by measuring k. 

A very recent experiment on driven domain walls [ j57| 
shows that the exponent k can be measured very accu- 
rately. This should make it possible to distinguish exper- 
imentally between both phases in the random bond case. 
For periodic media, an experimental realization is given 
by charge density waves. Recently measured TV curves of 
the conductor o£TaS3 at temperatures below IK can be 
fitted by Eq. fll3^ ) with k = 1.5 - 2 [Q. The experimen- 
tally observed tende ncy t o larger k for purer crystals is in 
agreement with Eq. (133) for D = 3. 



8 Discussion and Conclusions 

In this paper we have studied the static behavior of ori- 
ented elastic objects subjected to both a random potential 
and a periodic potential from an underlying crystal lat- 
tice. Scaling arguments and a functional RG calculation 
in D = 4 — e dimensions led us to the conclusion that the 
competition between both pinning effects induces a con- 
tinuous roughening transition between a flat phase with 
finite displacements and a rough state with diverging dis- 
placements. While based on the same physical mechanism, 
the description of the transition turned out to be different 
for periodic elastic media and single elastic manifolds. 

For periodic elastic media such as charge density waves 
or flux line lattices in high-T c superconductors, the transi- 
tion exists only above some critical ratio p c — l/ao of the 
lattice constants of the medium and the crystal potential. 
In D = 3 we obtain p c = 6/n from the e-expansion. Below 
Pc a weak random potential is irrelevant and the medium 
remains flat. On the flat side, this transition is described 
by three independent critical exponents which obey the 
scaling relations for random field systems with a discrete 
symmetry of the order parameter. Although the situation 
is similar to the thermal roughening transition in two di- 
mensions, where also logarithmic roughness exists in the 
rough phase, in the disorder driven case the existence of 
one universal stable fixed point describing the rough phase 
leads to a different RG flow. In contrast, for the thermal 
case the temperature axis represents a line of fixed points. 
Notice that in the disorder dominated system studied here, 
the important logarithmic roughness sets in only beyond 
the Fukuyama-Lee length La, whereas thermal roughness 
in two dimensions appears on all length scales. 
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In the case of elastic manifolds, the interplay between 
random potential and lattice pinning leads also to a con- 
tinuous transition for 2 < D < 4. We have obtained the 
exponent vn = for the divergence of the correlation 

length, which has to lowest order in e the same form as 
that of the thermal roughening transition if e = 2 — D. 
For both random field and random bond systems, elastic 
manifolds show a superuniversal logarithmic roughness at 
the transition. 

For D — 2 the interface is always asymptotically rough. 
At the upper critical dimension D = 4 an arbitrarily 
small lattice pinning produces a flat interface due to sub- 
logarithmic roughness in the absence of lattice effects. A 
diverging correlation length £|| appears on both sides of 
the transition. On the rough side £|| sets the length scale 
for the crossover from logarithmic to algebraic roughness. 

Recently, the exact ground state of interfaces in the 
three dimensional random bond, cubic lattice Ising model 
has been studied numerically |59[ . For interfaces oriented 
along the {100} direction, a roughening transition at finite 
disorder strength from an almost flat to a rough phase has 
been observed whereas interfaces along the {111} direc- 
tion remain always rough. Although our scaling arguments 
yield even for interfaces along the {100} direction no tran- 
sition, one has to take into account the finite system sizes 
used in the numerical calculation. Indeed, we expect a flat 
behavior of the interface width up to the scale which 
is exponentially large for weak disorder in D = 2, hence 
probably leading to a flat 'phase' at sufficiently weak dis- 
order in finite systems. 

More recently, also the replica Hamiltonian for periodic 
media with Ro((f>) ~ cos(</>), N = 1, has been studied us- 
ing the Gaussian variational method for p — 1,2<D<4 
|]60| and for p = 1/2, D = 3 Jll| . Both analysis yield a self- 
consistency equation for the effective mass of the field <^(x) 
in terms of the disorder strength. Solutions of this equa- 
tion exist only for v$ above a particular value vq, c , which 
is then identified with the phase boundary. At this phase 
boundary the effective mass remains finite corresponding 
to a first order transition. Moreover, the threshold «o,c is 
non- vanishing in spite of p < p c . As we have shown by scal- 
ing arguments and the RG calculation, these results can- 
not be correct in the case of weak disorder. To understand 
the origin of these inconsistent results, we note that the 
derivation of the implicit equation is based on the assump- 
tion of a replica symmetric solution for the rough phase, 
hence corresponding to the perturbative result of Larkin. 
In the rough phase, for the correct logarithmic roughness 
on scales beyond the Larkin length L^, replica symmetry 
breaking is needed within the variational approach. Since 
this logarithmic scaling on asymptotic scales is essential 
for the existence of a critical p c , the replica symmetric 
determination of the effective mass cannot reproduce our 
RG result. In addition, the self-consistency equation can- 
not expected to be valid in the vicinity of a second order 
transition since it does not capture the physics on scales 
^> La correctly. On the other hand, it should be noted 
that the variational approach, in contrast to the RG, does 
not depend on the assumption of weak pinning, hence de- 



scribing the system perhaps correctly in the case of strong 
pinning. However, the variational approach remains an un- 
controlled approximation. 

For p < p c , the correlation length £|| has to diverge for 
^0 - * vo.c = 0. Now we address the question of the cor- 
responding exponent z^|| . For this transition, there exists 
no unstable fixed point of the RG equations. As a first 
step towards an answer of this problem, one can calculate 
self-consistently an effective mass. To extend the valid- 
ity range of the self-consistency equation derived in Refs. 
[ p0]Jl3| to weak disorder, we have taken into account the 
correct asymptotic roughness by using the RG result for 
the case where vq = . Then the solution for the effective 
mass depends on p. For p < p c the mass remains always 
finite for vq > and vanishes continuously as vq = is 
approached, leading to the exponent v\\ = 1/2(1 — p 2 /p 2 c ) 
fell . If p > p c , the mass remains finite until vq. c > 0, 
where a solutions ceases to exist. In the latter case our 
RG treatment applies. We believe that the second order 
nature of the transition for p < p c is correctly given by our 
self-consistency condition, but the exponent estimated 
from this condition should not be taken seriously. 

The preceding discussion can be summarized in a spec- 
ulative phase diagram, which includes also the strong pin- 
ning case, see Fig. |. If the pinning potentials are strong 



enough such that the characteristic scales w 



-1/2 



and 



L, 



-l/(4-D) 



become of the order of max(ao,5o)j & 



first order transition may appear. 



(a) 



b. 




Fig. 9. Speculative phase diagrams including also the case of 



strong pinning where w 



-1/2 



and La 



A d - 1 become 



of the order of lo =max(ao,£o) for (a) p < p c and (b) p > p c . 
In the shaded region the RG is expected to be valid. Dashed 
lines represent possible locations of the roughening transition 
which may be of first order outside the validity region of the 
RG. 



In our study of periodic media we have not attempted 
to include the effect of topological defects. These can be 
considered if we treat </>(x) as a multi-valued field which 
may jump by ±27r at surfaces which are bounded by vor- 
tex lines. Including these defects, qualitatively we expect 
the following picture. For S = = vq it has been argued 
recently that for weak enough disorder strength Aq < Ajj , 
the system is stable with respect to the formation of vor- 
tices fj2| , |63 64 . However, vortex lines will proliferate for 
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Aq > Ar>. At present it is not clear whether the corre- 
sponding transition is continuous or first order. For 5 = 
but no > we expect that this transition extends to 
a line Ad(vq) until vq reaches a critical value vr> with 
Ad{vd) = A c (vd) (see Fig. [z]). For larger v the tran- 
sition is probably in the universality class of the p-state 
clock model in a random field, which has an upper critical 
dimension d c = 6. For p — 2 (Ising like case) the transi- 
tion is of second order |^5[ . A non-zero value of 6 will in 
general increase the size of the incommensurate phase, as 
schematically sketched in Fig. ||. 
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